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Abstract 


A model is developed to predict the shape of the distribution of 
granular agricultural particles applied by aircraft. The particle is 
assumed to have a random size and shape and the model includes the effect 
of air resistance, distributor geometry and aircraft wake. General 
requirements for the maintenance of similarity of the distribution for 
scale model tests are derived and are addressed to the problem of a non- 
general drag law. It is shown that if the mean and variance of the 
particle diameter and density are scaled according to the scaling laws 
governing the system, the shape of the distribution will be preserved. 
Distributions are calculated numerically and show the effect of a random 
initial lateral position, particle size and drag coefficient. A listing 
of the computer code is included. 
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particle drag coefficient 
aircraft lift coefficient 
particle diameter (m) 

2 

acceleration due to gravity (m/sec ) 

non-dimensional acceleration due to gravity 

constant in drag relation 

constant in initial velocity relation 

particle mass (kg) 

time (sec) 

aircraft speed (m/sec) 
non-dimensional x-velocity of particle 
x-velocity of air (m/sec) 
particle x-velocity (m/sec) 
non-dimensional y-velocity of particle 
y-velocity of air (m/sec) 
initial y-velocity of particle 
particle y-velocity (m/sec) 

magnitude of relative particle-air velocity (m/sec) 

non-dimensional z-velocity of particle 

z- velocity of air (m/sec) 

initial z-velocity of particle 

particle z-velocity (m/sec) 

longitudinal particle coordinate 

lateral particle coordinate 


y 



ii 


lateral ground intersection point 
initial lateral coordinate of particle 

z vertical particle coordinate 

initial vertical coordinate of particle 
a modified ballistic parameter 

3 ballistic parameter 

5 non-dimensional particle diameter 

Tj non-dimensional y-velocity of air 

p, mean (expected) value 

2 

V kinematic viscosity of air (m /sec) 

Q non-dimensional z-velocity of air 

p non-dimensional density of particle 

3 

p air density (gm/cm ) 

3 

p particle density (gm/cm ) 

P 

a square root of variance (standard deviation) 
T non-dimensional time 
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I. INTRODUCTION 


Aerial application of material has played an Increasingly significant 
role in United States agriculture, forestry, and other related industries 
since the Second World War. One of the major problems associated with the 
use of aircraft in this field has been the difficulty of obtaining desired 
distribution patterns with present equipment. This problem has been 
addressed by several groups in the academic and research fields, notably 
Henry (1962) and Yates, et al. (1970). However, most of this work has 
been of a trial and error nature with only sketchy guidelines being offered 
by theoretical considerations, mostly with regard to particle trajectories. 

It would appear, then, that an investigation of the factors which 
go to make up the structure of the distribution and which directly affect 
its shape would be of some use in designing better aircraft/distributor 
systems. Accordingly, this paper presents a model of the distribution of 
dry agricultural material which is based on the application of basic 
probability theory. The model takes into account the random nature of the 
size and shape of the particles, air resistance, distributor geometry and 
the effect of aircraft wake. Some results, based on data obtained for 
three grains (wheat, corn and oats), are presented in order to demonstrate 
the applicability and to show the effect of certain parameters on the 
distribution. In addition, a discussion of the effects of scaling on the 
distribution is presented and requirements for maintaining similarity in 
the distribution are discussed. 
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II. PARTICLE CHARACTERISTICS 

One of the first requirements in the development of a method for 
predicting the distributions of agricultural particles is a knowledge of 
the physical and aerodynamic characteristics of these particles. A 
literature search revealed that the necessary data were available for 
some particles but would have to be collected from numerous individual 
sources and investigations. An attempt was made to gather as much 
Information as possible on these characteristics from the literature. 

No attempt was made to make additional measurements . 

It was found that fairly complete characteristics could be assembled 
for several of the grains (corn, oats and wheat), while the same infor- 
mation was almost totally lacking for all of the solid fertilizers. The 
data presented here were collected from references (3-9), and were reduced 
via the method outlined below. 

The delineation of physical characteristics which best suited the 
point of view of this study was to assume the particle to have three 
fundamental characteristics which define its physical properties. They 
are; density, size (which together determine the weight) and shape (which 
determines the drag coefficient relation) . A sample collection of one 
type of particle will contain specimens possessing a wide range of each 
of these characteristics . That is, the characteristics will possess an 
apparent random nature, even though each individual specimen will have a 
fixed density and size, and will possess a definite relation between its 
drag coefficient and Reynolds number. The form of this relation is 
dependent on the shape of the particle (just as- a streamlined body possesses 
a different drag coefficient from a flat plate or a sphere) and thus 
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allowing the drag coefficient to be a random variable will reflect the 
random nature of the particle shape. 

Here the drag coefficient is defined in the conventional way 



where is the magnitude of the particle velocity relative to the air 
around it and S is a characteristic area associated with the particle, 
commonly taken to be its frontal area. As stated before, the drag co- 
efficient is a function of the Reynolds number of the particle, given as 



where d is a characteristic linear dimension of the particle. Above, p 
and V are the density and kinematic viscosity of the air, respectively. 

With this approach, then, the actual particle shape and its random 
nature is important only to the extent that it affects the drag coefficient 
relation and hence is completely independent of the size of the particle. 
Therefore, any convenient and consistent method for sizing the particles 
is acceptable. The approach taken here is one proposed by Garrett and 
Brooker (1965). The particle size is determined by considering a sphere 
which has the same volume as the particle in question and using the 
corresponding diameter of that sphere as the characteristic linear dimension. 
The characteristic area then is simply the cross-sectional area of the same 
sphere. The drag coefficient, measured by any of several methods described 
in the literature (the most common being a measurement of the terminal 
velocity of the particle), then corresponds to a Reynolds number calculated 
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using the diameter of an equal volume sphere. Thus, a drag 
relation can be constructed by a statistical analysis of the 
data. 

A least squares method was used to fit a general relation 
of the form 

Cj) = A Re + B/Re + C (1) 

to the data and the results are shown in Figures 1-3. This 
relation was taken to be the mean, or expected value, for the 
drag relation. A log-normal distribution function was used to 
describe the random variation of Cq about this mean. This dis- 
tribution has two primary properties : First , if a random var- 

iable X is lognormal, then In X is normal, or gaussian. 

Second, the function goes to zero as X approaches zero and is 
undefined for X less than zero. We should not expect any 
particles to have drag coefficients less than zero. 

In addition, this distribution, like the gaussian dis- 
tribution, is completely characterized by its mean value and its 
variance. The small amount and large scatter of the data for 
the variance makes any attempt at curve fitting questionable, 
however a least squares fit of a cubic 

Or = D Re? + E Re^ + F Re + G (2) 

was made. The results are presented in Figures 4-6. 

The particle size, as determined by the diameter of an 
equal volume sphere, was also held to be described by a lognormal 
distribution. However, the analysis was of a more straightforward 


nature and the results are 
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Reynolds Number 

Figure 5. Standard Deviation for 



standard Deviation of C 



Figure 6, Standard Deviation for Oats 
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tabulated in Table 1. It was found that the density of these particles 
did not vary significantly from specimen to specimen (probably for 
biological reasons) and so it is safe to assume the density to be a 
deterministic quantity. The values used are also presented in Table 1. 

In the calculation of the standard deviation of the drag coefficient, 
each point on the graph represents at least six data points, as rt was 
necessary to group neighboring points on the drag curve in order to obtain 
a standard deviation for a particular range of Reynolds number. 
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Table 1 


Physical Characteristics 



Diameter 

(microns) 



Mean 

Standard 

Deviation 

3 

Density (g/cm ) 

Oats 

3594 

426 

1.369 

Corn 

8133 

338 

1.234 


Wheat 


3573 


189 


1.320 
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III. TRAJECTORY EQUATIONS 


Each particle, as it is ejected from the aircraft, becomes a free 
projectile and is subject to Newton's second law of motion. The equation 
of motion, in vector form, describing the trajectory is 


m + D+ L+ mg = 0 


( 3 ) 


where L refers to the lift force generated by the particle. Since the 
particles studied here are irregular in nature (in particular, oats have a 
flattened profile), they will, in general, produce a force perpendicular 
to the direction of motion, as well as one opposite the direction of motion. 
These forces are termed lift and drag, respectively. Though, as shown in 
the preceding section, there is information regarding the drag experienced 
by these particles, predicting and representing the lift forces is a more 
complex matter. It is felt that, by ignoring the lift force, the sub- 
sequent simplification will be of such a degree as to offset the relatively 
small loss of accuracy. 

With this simplification, then, the single vector equation can be 

written more revealingly as three scalar equations, with the component of 

til 

the drag force in the i— direction being 


D. 
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Pi 


u. 


- U 


ia 



so that the drag acts in the opposite direction to the velocity vector 
' du 



+ D = 0 
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(A) 



m 




dt 


dw 


m -T“ + D - mg = 0 

utl Z 
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(5) 

( 6 ) 


It is apparent, however, from the work of Bragg (1977), that the 
motion of the particle remains predominantly in a plane perpendicular to 
the direction of motion of the aircraft (the x-direction) and since the 
distribution is dependent only on the motion in this plane, a two 
dimensional model will be used for this study. 

Rewriting the y and z equations of motion with the three fundamental 
particle characteristics explicitly shown, gives, after dividing through 
by the particle mass 


dv _ 


dt 


^ Pr 


(v -V ) 

p 


!]_ l/ (v -V )^ + (w -w )^ 
D “ p a 'pa' 


= 0 


(7) 


dw 

_£ _ 

dt 






V (v )^ + (w -w = 0 . 
D * ' p a' 'pa' 


( 8 ) 


Since the wake system of the aircraft is dependent primarily on the 
aircraft geometry and velocity, the nondimens ionalization was carried out 
with respect to the aircraft flight speed U and the wxng semi- span, b. 

The equations then become 

S + + ('*-0^ - 0 (9) 

I? - s* + 1 -r S + (w-o^ = 0 (10) 

> 

where the nondimensional variables are defined in Table 2. 
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Table 2 

Nondimensionalization 


v,w 


V 

U 


w 

_E 

U 


n»c 
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Detailed numerical studies of the above system have been carried out 
by Reed (1953) and Bragg with respect to liquid droplets. Although the 
trajectories predicted in both cases seemed to indicate good agreement 
with experiment, these analyses could only yield distributions by the time 
consuming and costly method of a numerical compilation of trajectories 
calculated using a matrix of initial conditions and particle parameters. 

It seemed desirable to construct an approximate solution to the above 
system which would offer a simple method of constructing the distributions 
while preserving some of the accuracy of the previous analyses. 

Returning to equations (9) and (10) , one of the most striking aspects 
revealed there is the fact that the particle-dependent parameters can be 
grouped into a single variable, and it is this variable alone which affects 
the solution of the system. This parameter bears a close resemblance to 
the ballistic coefficient, which is used in the study of projectile motion, 
and so it can be referred to as the ballistic parameter, p, where 


o = 3 fa fo 
P - 4 pp 6 


Because of the presence of the drag coefficient, the ballistic parameter 
is a function of the particle Reynolds number. However, extensive 
simplification of the entire system can be achieved by observing that a 
simple drag relation of the form 

Cjj = K/Re (11) 

offers nearly as good a fit to the data as equation (1) . 

Substitution and further manipulation with this result yields the 


following system: 
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+ a(v-Tl) = 0 (12) 

^ - g* + a(w-Q = 0 <13) 

where a modified form of the ballistic parameter, a> has been used 


3 iL fs ^ 


(14) 


Although it appears that much simplification has been achieved, this 
system still requires numerical solution due to the complex nature of the 
dependence of T) and Q on position. 

The range of values assumed by a for the grains studied in this 
report necessitated the use of two approximate solutions to this system. 
The first is to consider a to be vanishingly small. This corresponds to 
the case where the particle is quite large and dense, or the medium 
through which it travels is quite ratified and the drag coefficient 
low. This results in the classical case of a body in free fall in 
a constant gravitational field. The solution for the trajectory is 
simply 


y ~ y + v T 
^ o o 

(15) 

2 

Z = z +WT-g*— . 

o o 2 

(16) 


The lateral position at which the particle strikes the ground can be 
calculated by letting z = 0 and solving equation (16) for t, and sub- 
stituting into equation (15). The result is 

y^ a y + — X (w + Jw^ + zz g*) . 

•'G -^o g* ^ o Vo 0 ° ' 


( 17 ) 
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Comparison of this result with a numerical integration of equation (12-13) 
showed excellent agreement for values of a up through about .20 
corresponding to a diameter of about 3500 microns, as shown in Figure 7. 

The second approximation was arrived at by considering initially 
the case of a particle falling through still air, that is, neglecting the 
aircraft wake. A study of this system revealed that the smaller particles 
(larger a) tended to be carried further outboard' and also have a 
significantly lower flight time in the numerical computation than in this 
approximation. It is apparent that the effect of the wake behind the air- 
craft is to draw the particles downward and outward from the aircraft 
fuselage. This effect is primarily due to the pair of vortices whxch are 

shed from each wing tip as a result of the lifting action of the wing and 

1 

which cause the air in the vicinity of the vortices to rotate in the 
manner shown in Figure 8. 

In accordance with these observations, a modification to the still- 
air model was made which extended the range of validity to cover the grains 
being studied. This modification consists of adding a constant air 
velocity whose magnitude and direction is consistent with those found in 
an aircraft wake and which result in the best agreement with the numerical 
calculation. The adjusted equations then become 

+ ccv - all = 0 (18) 

^-g* + ow-aC=0 (19) 

where T| and Q are constants. The solution for the trajectory is simply 

V -T| 

y “ + T]t + (l-e’'^'^) (20) 
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y (semispans) 
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Figure 7. Approximate Solutions 




Figure 8. Wake Structure 
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( 21 ) 


Sample comparisons between the trajectories calculated with equations 
(15-16) and (20-21) and those arrived at numerically are shown in 
Figures 9-14. Since the cases with larger a had longer flight times, an 
expression for the lateral position of the ground Intersection can be 
derived by assuming 



« 0 


resulting in the expression 


= y + 
o a 




z a + w 

+ j (-2-- 2) 

a Si . r 

a 


( 22 ) 


Comparison of this result with the numerical calculation is also 
shown in Figure 7. It will be seen later that this extensive approximation 
is necessary in order to conduct the probability analysis . 
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IV, THE DISTRIBUTION MODEL 


As stated earlier, this model makes use of results from basic 
probability theory. 

In general, the ground Intersection point, y , of a particle which Is 
ejected from an aircraft Is a function of the Initial conditions of its 
trajectory. Its particle parameters, the aircraft geometry and flight 
conditions. For the system defined here, with the appropriate non- 
dimenslonallzatlons, this can be written as 




( 23 ) 


where C , the aircraft lift coefficient, Is a measure of the magnitude of 

Li 

the wake velocities induced by the tip vortices. In the sense that these 
variables can assume a range of possible values for any given type of 
particle and aircraft/distributor system, they can be thought of as random 
variables. As such, they will possess a joint probability density function 
(PDF) characterizing their random nature of the form 



,Z ,V ,W ,A, 
o o o 


L 


This is the most general representation of the probabilistic system in 

t 

that the individual PDF's can be derived from it by first fixing the values 
of the other variables, then multiplying by the probability that each will 
assume that fixed value, and finally integrating over all possible values 
of the fixed variables, so Chat, for instance. 



O 3 V W OCt o' o’ o o' ’ L 
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W « VT 

o o o L 


?< P„ „ „ . - (2 ,v ,w -CL.C, )d2 dv dw dadC^ 

Z »V ,W ,A,C. o’ o’ o’ L o -o 0 L 
o o o L 


or more simply 


% " I ••■ L ,Z ,W ,A,C 


2 v^w^aC^ o’ o’ "o’ *L 


X dz^dv^dw^dodC^ 


Purthar, since it is known that if a function of & random variable defines 
a one to one laapplag onto another random variable, the PDF for the second 
variable is defined as 


Py(y) “ 


l|! 


so that by again fixing all the variables in equation (23), an expression 
for the PDF of y„ can be derived. Writing first 



then multipLying through and integrating as before, the result is 


y (Fq) = J ... J 2 ,v ,W ,A,Ci 

G ,w ,a,C, o’ o’ o’ o’ 

O O D id 


^ G 


( 24 ) 
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Considering that many particles go to make up a single deposition, the 
probability density function will give a good representation of the actual 
deposition shape. Note also that the same result can be arrived at by 
fixing any combination of the variables and integrating over those. In 
addition, since it can be assumed that the particle ballistic parameter 
is independent of the initial conditions, that is, the particle character- 
istics do not depend on the distributor geometry, the joint PDF" can be 
split and equation (24) can be rewritten as 



%) 




- f, ,Z ,V ,W ,C ^^o’*o’^o’”o’^L^ 



(25) 


This, then, is the general equation describing the ground distribution. 
The necessity of an analytical expression for y in terms of the other 
parameters is apparent due to the presence of the term dy /dy . 

O Vy 
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V. SOME DISTRIBUTIONS 


In the interest of demonstrating the applicability of the model 

derived in this study, only two parameters were considered to have a 

random nature, y^ and a» The other parameters were represented as 

functions of these, or as constants. Extension to the more general case 

is straightforward and merely requires additional integration. 

Specifically, z and C,. were held as constants, and v and w are 
O o o 

represented by 

\ = Vo 


The resulting simplified expression for the PDF of y is 

\ G 


G ^ o G 


(26) 


where 


dy 


1 + — (w + /w^ + 2z g*) 
— o ’ o o 


Idy 


1 +^ 


0 < a < . 20 


a > .20 


(27) 


a 


and the PDF for a can be derived as follows. Equation (14) can be 
rewritten as 


In a “ In C + In K - 2 In 6 . (28) 

A basic result of probability theory is that the sum of any number of 
gausslan random variables is also gaussian. It can then be concluded 
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that, since In K and In 5 are gaussian. In a is gauss ian. The effect of 
the constant is merely to shift the mean. Finally, since In a is 
gaussian, a is lognormal. Its mean and variance are computed as follows. 
Let 

Z = In a 
Y = In K 
X = In 6 
C* = In C . 

Then 

Z = C* + Y - 2X 

and the expected value of Z, or the mean, written as E[z], is simply 

E[z] = C* + E[Y] - 2E[Z] , 
a more common notation for the mean is |j,, so that 




C* + Uy - . 


2 

The variance, written as eCCZ-ij^,) ] is simply 


(29) 


or 


E[(Z^-2Zp^+p^^)] = E[Z^] - 2E[z]p^ + 


Substituting for Z and expanding as above 


0^^ = E[Y^] + 4e[X^] + C*^ + 2c*E[Y] - 4C*E[X] - 4E[Y]E[X] - . 
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Further substitutloa gives 


+ (C* - 2yjj + 


or 


2 2^,2 

°z =°Y 


(30) 


The geaeral formula for the lognormal PDF is 


P^(t) = 

CTt 


t - p.) ^ t > 0 

2a 


(31) 


vhere the values for a and pi are those used above. The actual expected 
value of t is given by 

00 

E[t] = J t • P^(t)dt 
0 

which results in 

^ . E[t] - .1* + T- , , (32) 

2 

and the variance of t, (j^ is given by 


CO 

= E[(t-p,^)^] = J t%^(t)dt - p^^ 
0 


and the result is 






(33) 


Equation (32) and (33) furnish relations between the actual mean and 
standard deviations of a random variable t (p.^, and ct^) and those values 
used in the lognormal PDF (p, and a). With these results, then, the PDF 


for a can be written as 



ORIGINAL PAGE IS 
OF POOR QUALITY 


34 


p (( 3 ^) = — exp [ ^ ~ ] (34) 

^ aa ^ 2a 

where 

li = In C + 

2 2,2 
CT = + 4o2 

and 

a/ . InCl + (^fy 
1 

^ - In ^ - -j- 

, 

= In - — • 

A computer program was written to numerically integrate equation (26), 
applying the above results. It is called DEP and a listing is provided 
in the Appendix. The variation of y^ is provided in the form of FUNCTION 
subroutines PYOF and PYOG which can be easily modified to account for any 
type of variation of y^. For the case in which y^ is a constant, that is, 
there is a single orifice or nozzle, the resulting distribution for wheat 
is shown in Figure 15. 

For the case in ^ich the PDF of y^ is constant across a given width, 
results are shown in Figures 16-20 in which the type of material, then the 
width of the distributor is varied. 



(yp) (l/semispans) 



Figure 15 


Wheat, .05 jf ^ *06 
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(y ) (1/semispans) 




(y„) (1/semispans) 



Figure 17 


Wheat, 0 < y <0.2 
> — •'o — 
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VI. SCALING CONSIDERATrONS 

In order to facilitate experimental verification of these results and 

also to aid in obtaining viable results from test programs where scale 

1 

model aircraft are used, it would be useful to investigate the effect of 
scaling the entire system on the probabilistic aspects of the problem. 
Specifically, it would be desirable to be able to predict probability 
density functions for the scaled values of the parameters in a given test 
which would result in distributions which are geometrically similar to 
those in a corresponding full scale system or a test of a different scale. 

The problem of obtaining proper scaling for a deterministic system 
has been addressed by Ormsbee and Bragg (1978) and the derivation will not 
be presented here. Instead, a heuristic argument will be presented in 
order to gain insight into the nature of the probability problem, and then 
a formal approach will be undertaken. 

Consider a collection of particles which are to be ejected from an 
aircraft. If two individual specimens are singled out, one possessing the 
mean diameter of the collection and the other possessing a diameter which 
is smaller by one standard deviation, and the trajectories of these 
particles are mapped, they would appear similar to those in Figure 21(a). 
In this argument, all other factors (i.e., the drag relation, the initial 
conditions and the density) are held to be deterministic. In order to 
obtain trajectories which are geometrically similar to each of these from 
a scale model test, it would be necessary to obtain two specimens from a 
collection of a different type of particle such that the requirements 
for scaling are met. These particles would describe trajectories similar 
to those shown in Figure 21(b). 
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It seems to follow that If the second collection of particles 
possesses mean and standard deviation values corresponding to the diameters 
of the specimens chosen, scaling of the distribution will have been 
achieved if the nondimens ional lateral ground positions indicated xn 
Figure 21 are identical. This presupposes that the shape of the prob- 
ability density function has not changed so that it is sufficient only to 
insure that its width (a measure of which is the variance) has remained 
the same. This can be shown to be true if the scaling laws are observed. 

More formally, proper scaling of the distribution will obtain if the 

probability density function for the modified ballistic parameter, a, 

remains identical. The only other possible inputs to the system are the 

geometry and the wake system and these have been accounted for by the 

deterministic scaling laws. A probability density function is completely 

th 

determined by the values of its moments, the i — moment about the mean 
being given by 

^^i P3j(x)dx . 

X 

Since there are an infinite number of moments associated with a PDF, in 
general it is not possible to insure complete similarity. However, for 
the special case of the gaussian distribution (from which the lognormal is 
derived), all of the moments are not independent and, in fact, all the 
higher order moments can be written in terms of the fxrst two: the mean 

and the variance. Thus the problem of distribution scaling is reduced 
to the question of maintaining the mean and standard deviation of the 
ballistic parameter. 
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In order to retain as much generality as possible, the earlier 
assumption that the density is deterministic will be removed and it will 
be assumed to have a lognormal PDF. Equation (28) can be rewritten to 


include this change as 

In a = In C ' + In K - In p - 2 In 6 (35) 

where p is a nondimensional density obtained by dividing the particle 
density by the density of air. Letting 

r = In p 

and proceeding as before, the mean and variance of a can be determined, 
resulting in 

= e (36) 

where 

H = In C' + py - lij. “ 2y,^ 

2 2 ^ 2^,2 

a - a + cf + 4a 

y r X 

and 

2 1 n /^K.2 t 

°y “ ’ 

a ln[l + (--) ] 

* ln[l + 



45 


2 

2 

2 

2 

lij. - In p,p - 2 

Thus scaling will be achieved If the values of the mean and variance 
as given by equations (36) and (37) are the same for the scaled system as 
for the full scale. 

In the case of a general drag law, the deterministic values for the 
particle diameter and density necessary to insure scaling, if the aircraft 
is reduced in size by a factor of s, are as follows: 


d = d- Ts 
s f ^ 



The corresponding nondimens ional variables will have the following values: 


6 


s 



3 

2 


Ps = 



For this case it is necessary for the scaled particle to retain the same 
Reynolds number as the full scale and for the drag law to remain identical. 
The values required for the mean and variance of the scaled diameter and 
density which will insure scaling can be derived in the manner below. 
However, since both transformations are of the form 
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X = X- s 
s f 

it will suffice to carry through the derivation for just the diameter. 
First, observe that 

3 

In 6 = In 5 . + -T In s . 

S X ^ 


Since In 6- is gauss ian, this demonstrates that In 6 must be gaussian 
X s 

and, consequently, 5 must be lognormal. Letting 

s 


and rewriting 


Z = la 6^ 

X ■= In 6- 

C = In s 
Z = X + C . 


The mean and variance of Z can be computed as before, giving 


= tS, + C 


2 2 

a ~ <j 

Z X 


and, since 6 is lognormal 
s 




- e 




+ C + 


X 


= e n- 
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or 



Likewise 


or 



1 ) 



1 ) 




The derivation for the density proceeds identically, and results in 


-3/2 


Pf 


-3/2 
a. = 3 a 


Pf 


Substitution back into equations (36) and (37) will verify the 
scaling since 


C ' 
s 



a 


s 


2 



P-s 


In C' 



In 3 + u “ 


+ I 1» s 



- 3 In s 
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If the analysis is restricted to the case where a simpler drag 
relation of the form 

C 

Re 

is used (as in this study), a different result is obtained. As discussed 
by Ormsbee and Bragg, an approximation such as this frees the investigator 
from the separate requirements on density, diameter and drag relation and 
allows him instead to consider only their combination in the ballistic 
parameter. Similarly, in the probabilistic situation, it is not necessary 
to meet the requirements in the analysis immediately preceding, but only 
the more general ones defined by equations (36-37). Thus, as long as the 
PDF's for the density, diameter and drag relation are lognormal, any 
combination of means and variances for these parameters which satisfy 
equations (36-37) could be used in a scale test and still preserve the 
shape of the' distribution. 
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VII. CONCLUSIONS 


For granular particles which are larger than approximately 2500 
microns in diameter, considerable simplification of the trajectory equations 
can be achieved due to the limited influence of the aircraft wake. The 
use of an approximate drag relation, which allowed further simplification, 
does not restrict the applicability as much as would at first appear, 
since numerical calculations show that the particle Reynolds number is 
maintained near the value associated with its terminal velocity for the 
entire time span of the trajectory, except for an initial acceleration 
which is of short duration. Thus the drag relation used need only 
approximate a small part of the drag curve in the range of the terminal 
Reynolds number. 

In the case of smaller particles (below 2500 microns) the probability 
analysis could still be used if an empirical relation of the form 


Fg = f(7o.a) 


is used instead of the algebraic expressions derived here. 

dy 

necessary to compute the Jacobian of the transformation, 

G 

able to write the above expression in a form 


It is only 
, and to be 


= g(7p»a) 
o G 

so that the distribution of y can be calculated. 

The usefulness of this analysis can be demonstrated by noting that a 
distribution can be calculated by a single integration, whereas a 
compilation of trajectories would require one integration for each 
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trajectory. Moreover, if more than one initial condition is allowed to 
vary, the cost of producing a distribution by this method will increase 
by a factor equal to the number of initial conditions varied. On the 
other hand, the cost associated with a compilation of trajectories will 
increase by the power of the number used — a geometric, rather than an 
arithmetic, increase. 

It was shown that the shape of the distribution for a scale model 
test will be similar to that of a full scale test if, in the most general 
case, the mean values and variances of the particle parameters are scaled 
just as the deterministic values would be. Moreover, it was shown that 
if a drag relation of the form 

C 

D Re 

is used instead of a more general one, similarity could be maintained in 
any number of ways simply by observing a restriction on a combination of 
the means and variances of the particle parameters, rather than requiring 
separate concurrence with the scaling laws. This simplification offers a 
greater freedom in the choice of particles to be used in tests verifying 
these results and in any further work in which scale models are used. In 
view of the fact that particles meeting the requirements of the more 
general set of scaling laws have proved to be rather difficult to locate 
and/or implement, the simplified procedure for scaling offers a versatility 
which, hopefully, will more than offset any inaccuracies involved. 
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APPEMDIX 

The computer code DEPOSIT was written to calculate the distributions 
by numerical integration using a PORTUOI catalog subroutine entitled 
GQU3Z. The input, where dimensions are required, can be done in any 
consistent system of units, due to the fact that all values are non- 
dimensionalized within the program, except the particle diameter, which 
must be input as the nondime ns ional 6, where 



A sample printout is shown, depicting the form of the output. By 
manipulating the variables RES and SC, the code can deliver any resolution 
of the distribution and can extend the range of calculation to cover any 
length of the lateral ground coordinate, y„. Various types of distributor 
geometries can be analyzed by modification of the FUNCTION subprograms 
PYOF and PYOG. PYOF covers the approximation that 

ttT » 1 


while PYOF assumes that 


a « 1 . 
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PROGRAM DEPOSIT (INPUT, OUTPUT) 

REAL NU,MD,MK,MA,MD1,MK1,MA1,KV 
EXTERNAL AUX 

COMMON/COM1/ YG,MA1,SA1,PI,MA 

C0tW0N/C0M2/ ETA , XI , KV , WO , ZO , G , YI , YF 

READ2,MARK 

PRINT4 ,MARK 

PRINTS 

READS, U,B,ZO,NU,RA 
PRINTS ,U,B,Z0,NU,RA 
READS, RP,MD,MK,SD,SK 
PRINT? , RP , MD , MK , SD , SK 
READS, YI,YF,KV, WO SC 
PRINTS , YI , YF , KV , WO , SC 
READ3,UP,H,RSS,ETA,XI 
PRINT80 , UP , H , RES , ETA , XI 
PRINTS 

PI=S. 14159265 
G=9.808«B/(U«U) 

C=0 . 75*NU»RA/ (U*B»RP ) 

SD 1 =SQRT ( ALOG (1 . 0+ ( SD /MD ) «» 2 ) ) 

SK1=SQRT(AL0G'( 1 .0+(SK/MK)»*2) ) 

MD 1 =AL0G (MD) -SD 1 »SD 1 /2 
MK1=AL0G(MK)-SK1«SK1/2 
PRINT10 ,C,SD1 ,SK1 ,MD1 ,MK1 
SA 1 =SQRT ( 4»SD 1 *SD 1 +SK 1 *SK 1 ) 

MA 1 = ALOG ( C ) +MK1 -2 . «MD 1 
MA=EXP(MA1+SA1»SA1/2 ) 

SA=SQRT (MA*MA* ( EXP ( S A 1 *SA 1 ) - 1 . ) ) 

PRINT11,SA1,MA1,SA,MA 

PRINT 12 

FYG=O.OSO 

NRES=IFIX(RES)+1 

D01 Isl,NRES 

YG=(I-1)»SC/RES 

CALL GQU3Z(0.0E0,DP,AUX,H,PYG) 

FYG=FYG+PYG»SC/RES 

PRINT13,I,XG,PYG,FYG 

1 CONTINUE 
STOP 

2 FORMAT(AIO) 

3 F0RMAT(5E10.4) 

4 FORMAT ("1", ///I OX, "GROUND DEPOSITION OF ",A10) 

5 FORMAT (//13X, "INPUT PARAMETERS") 

6 F0RMAT(/"U =" ,F5 . 1 ,9X, "B =",F5.2,9X,"Z0=",F5-3,9X,"NU=",E12.6,2X," 
&RA=",E12.6) 

7 F0RMAT(/"RP=",E12.6,2X,"MD=",E12.6,2X,"MK=",E12 6 2X,"SD=" ,E12.S,2 
&X,"SK='’,E12.6) 

8 F0RMAT(/"YI=" ,F5 . 3 , 9X , "YF=" , F5 • 3 ,9X , "KV= " , F7 . 5 , 7X , "WO=" ,F8 . 5 , 6X , "S 
&C=",F5.3) 

80 FORMAT (/"UP=",F6. 3, 8X,"H =",F7.3,7X,"RE=",F7.3,7X,"ET=",E12,6,2X," 
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&XI=",E12.6) 

9 format!//, 1 1X, '‘CALCULATED PARAMETERS") 

10 format!/" C =",E12'.6,2X,"SD1=",E12.6,2X,'‘SK1='*,E12.6,2X "MD1=",E12 
&.6,2X,‘‘MK1s",E12.6) 

11 F0RMAT!/"SA1 = '‘,E12.6,2X,‘‘MA1 = '',E12.6,2X,'‘SA r",E12.6,2X,'’MA =",E12 

&. 6 ) 

12 format! "1",///3X," I ",6X," IG ",9X,"PYG(YG)",8X,"FXG!YG)") 

13 F0RMAT!3X,I3,3X,E12.6,3X,E12.6,3X,E12.6) 

END 

SUBROUTINE AUX!A,PY) 

REAL MA1,MA,KV 

C0MM0N/C0M1 / YG , MA 1 , SA 1 , PI , MA 
C0MM0M/C0M2/ ETA,XI,KV,WO ,Z0 G,YI,YF 
IF!A.LT. 0.230) GOTO 1 
DYDG=1./!1.+KV/A) 

PY0=PY0F!YG,A) 

GO TO 2 

1 DYDGsl ./! 1 .+KV»!WO+SQRT(WO»WO+2»ZO»G))/G) 

PY0=PY0G!YG,A) 

2 PEX=!AL0G!A)-MA1 )»»2/!2*SA1»SA1 ) 

IF!PEX GT.600.) GOTO 21 
PA=EXP!-PEX)/!SA1»A»SQHT!2»PD) 

GOTO 22 

21 PAsO.OEO 

22 CONTINUE 
PYsDYDG»PY0»PA 
RETURN 

END 

FUNCTION PYOF!YG,A) 

REAL KV 

COMMON/COM2/ ETA, XI, KV, HO ZO,G,YI,YF 
Y0s!YG-ETA»!Z0»A+W0)/!A»!G/A-XD) )/(1 .+KV/A) 

IF!YO.LT.YI) GOTO 1 
IF!YO.GT,YF) GOTO 1 
PY0F=1./!YF-YI) 

GOTO S 

1 PYOF=O.OEO 

2 RETURN 
END 

FUNCTION PYOG!YG,A) 

REAL KV 

C0MM0N/C0M2/ ETA ,XI ,KV , WO , ZO , G , YI , YF 
YO=YG/ ! 1 . +KV» ! WO+SQRT ! W0*W0+2» ZO »G ) ) /G ) 

IF!Y0.LT.YI) GOTO 1 
IF!YO.GT.YF) GOTO 1 
PY0G=1./!YF-YI) 

GOTO 2 

1 PYOG=O.OEO 

2 RETURN 
END 



1 


u = 

RP= 

yi=( 

UP= 

c ; 

SAIi 

1 


GROUND DEPOSITION OF WHEAT 


INPUT PARAMETERS 


30.5 

B = 6.77 

Z0= .500 

NU= 

.146700E-04 

RA= .122500E-02 

.132040E+01 

MD= .528100E>.03 

MK= .127400E+04 

SD= 

279200E-04 

SK= . 2265 OOE+O 3 

i.OOO 

YF= 300 

KV= 1.00000 

W0= 

.10000 

SC=2.000 

2.000 

H = 48.000 

RE=100.000 

ET= 

.160000E-01 

XI=- I 8 OOOOE-OI 


CALCULATED PARAMETERS 

.494673E-10 SD1= .528319E-01 SK1= .176405E+00 MD1=-.754762E+01 MK1= .713436E+01 
. 20563 OE+OO MA1=-. 15001 1E+01 SA = .473571E-01 MA = .227873E+00 





I YQ 

1 0 . 

2 .200000E-01 

3 .400000E-01 

4 .600000B-01 

5 .800000E-01 

6 .100000E+00 

7 .120000E+00 

8 .140000E+00 

9 .160000E+00 


PYG(YG) 

.2907I3E+OO 

.29O7I3E+OO 

.290713E+00 

.600057E+00 

.603464E+00 

.603464E+00 

.603464E+00 

.603464E+00 

.603464E+00 


FYG(YG) 
.581426E-02 
.116285E-01 
.174428E-01 
.294439E-01 
.415132E-01 
•535825E-01 
. 65651 8E 01 
•777210E-01 
. 89790 3E-0J 


01 

ui 
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26 

.500000E+00 

.603464E+00 

. 294968E+00 

27 

.520000E+00 

.603464E+00 

•307037E+00 

28 

.540000E+00 

.603464E+00 

.319107E+00 

29 

.560000E+00 

.603464E+00 

.331176E+00 

30 

.580000E+00 

.603464E+00 

. 343245E+00 

31 

.600000E+00 

.603464E+00 . 

•355315E+00 

32 

.620000E+00 

.603464E+00 

. 367384E+QO 

33 

.640000E+00 

•603464E+00 

. 379453E+00 

34 

.660000E+00 

.603464E+00 

.391522E+00 

35 

.680000E<i-00 

.603464E+00 

.403592E+00 

36 

.700000E+00 

.603464E+00 

.415661E+00 

37 

.720000E+00 

.603464E+00 

.427730E+00 

38 

•740000E+00 

.603464E+00 

.439800E+00 

39 

.760000E+00 

.603464E+00 

.451869E+00 

40 

•780000E+00 

.603464E+00 

.463938E+00 

41 

.800000E+00 

.603464E+00 

.476007E+00 

42 

.820000E+00 

.603464E+00 

.488077E+00 

43 

.840000E+00 

.603464E+00 

.500146E+00 

44 

.860000E+00 

.603464E+00 

.512215E+00 

45 

•880000E+00 

.603463E+00 

.524285E+00 

46 

.900000E+00 

. 60346 1E+00 

•536354E+00 

47 

•920000E+00 

' .603457E+00 

■548423E+00 

48 

.940000E+00 

.603448E+00 

.560492E+00 

49 

.960000E+00 

.603429E+00 

•572560E+00 

50 

-980000E+00 

.603392E+00 

.584628E+00 

51 

.100000E+01 

.603329E+00 

.596695E+00 

52 

.102000E+01 

.603217E+00 

.608759E+00 

53 

.104000E+01 

.603032E+00 

.620820E+00 

54 

.106000E+01 

.602748E+00 

.632875E+00 

55 

.108000E+01 

.602319E+00 

. 64492 1E+00 

56 

.110000E+01 

.601672E+00 

.656955E+00 

57 

.112000E+01 

.600782E+00 

.668970E+00 

58 

.1140008+01 

. 59950 1E+00 

.680960E+00 

59 

.116000E+01 

.597838E+00 

.692917E+00 

60 

.118000E+01 

.595683E+00 

.704831E+00 

61 

•120000E+01 

.592904E+00 

•716689E+00 

62 

•122000E+01 

•589319E+00 

.728475E+00 

63 

.124000E+01 

, 58476 lE+00 

.740170E+00 

64 

.126000E+01 

. 57957 lE+00 

.751762E+00 

65 

.128000E+01 

•573285E+00 

.763227E+00 

66 

.130000E+01 

.566227E+00 

.774552E+00 

67 

•132000E+01 

.557848E+00 

.785709E+00 

68 

.134000E+01 

. 54821 3E+00 

.796673E+00 

69 

.136000E+01 

.536765E+00 

.807409E+00 

70 

.138000E+01 

.525223E+00 

.817913E+00 

71 

•140000E+01 

.512360E+00 

•828160E+00 

72 

.142000E+01 

.498569E+00 

.838132E+00 

73 

.144000E+01 

.482121E+00 

.847774E+00 

74 

.146000E+01 

.467171E+00 

.857117E+00 

75 

.148000E+01 

.451056E+00 

.866139E+00 

76 

.150000E+01 

.432445E+00 

.874787E+00 
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77 

.152000E+01 

.415287E+00 

.883093E+00 

78 

.154000E+01 

•396127E+00 

i8910l6E+00 

79 

.156000E+01 

.377620E+00 

.898568E+00 

80 

.158000E+01 

.358334E+00 

.905735E+00 

81 

.160000E+01 

.339880E+00 

.912532E+00 

82 

.162000E+01 

.318155E+00 

.918895E+00 

83 

.164000E+01 

•301217E+00 

.924920E+00 

84 

.166000S+01 

.290713E+00 

.930734E+00 

85 

.168000E+01 

.290713E+00 

•936548E+00 

86 

•170000E+01 

.29071 3E+00 

.942363E+00 

87 

.172000E+01 

.290713E+00 

.948177E+00 

88 

•174000E+01 

•290713E+00 

. 95399 1E+00 

89 

,176000E+01 

.290713E+00 

.959805E+00 

90 

.178000E+01 

.29071 3E+00 

.965620E+00 

91 

•180000E+01 

.29071 3E+00 

.971434E+00 

92 

.182000E+01 

. 29071 3E+00 

.977248E+00 

93 

.184000E+01 

.290713E+00 

.983062E+00 

94 

.186000E+01 

.290713E+00 

.988877E+00 

95 

.188000E+01 

.290713E+00 

.994691E+00 

96 

.190000E+01 

.29071 3E+00 

.100051E+01 

97 

.192000E+01 

0. 

.100051E+01 

98 

.194000E+01 

0. 

.100051E+01 

99 

•196000E+01 

0. 

.100051E+01 

100 

.198000E+01 

0. 

.100051E+01 

101 

.200000E+01 

0. 

.100051E+01 
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A Method for the Determination of the Matrix Y. K, Lin 
of Impulse Response Functions with Special 
Reference to Applications in Random 
Vibration Problems 


AFFDL-TR-66-80 . Wright 
Patterson AFB, 743-751 
(1966), 


Convergent Channel Shock Tube for Detonation A, J. Crooker 
Initiation Studies 


"Detonation and Initiation 
Behind an Accelerating Shock 
Wave" by R. A. Strehlow, 

A. J. Crooker, R. E. Cusey, 
Comb and Flame , 11:339-351 
(1967). 


A Comparison of Experimental and Theoretical 
Transverse Wave Spacings in Detonation 

A Simple Model for the Mechanism of Detonation 

Transverse Wave Structure in Detonations 

A Real Gas Analysis Using an Acoustic Model 
for the Transverse Wave Spacing in Detonations 

Shock Tube Studies in Exothermic Systems 


R, 

H. 

Watson 

See AAE 66-3 

J. 

R. 

Eyman 

See AAE 66-3 

R. 

A. 

Strehlow, 

11th. Symposium (Inter- 

R. 

Liaugminas, 

tional) on Combustion, 

R. 

H. 

Watson and 

Mono Book Corp. 

J. 

R. 

Eyman 

Baltimore, Md., (1967), 

R. 

E. 

Maurer 

AIAA Journal, 7; 323-328 




(1969). 

R. 

A. 

Strehlow 

Phys. Fluids, 12: 96-100 


(1969). 
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AAE 69-1 


Title Author 

Shock Tube Chemistry R. A. Strehlow 


Structural Failure Criteria for ■ A. R. Zak 
Solid Propellants Under Multiaxial 
Stresses 

Structural Analysis of Realistic 
Solid Propellant Materials 

Characteristics of Transverse Waves 
in Detonations of 
and CH^ - Oxygen Mixtures 

A Review of Shock Tube Chemistry 


On the Interpretation of Molecular 
Beam Data 

Detonative Mach Stems I 


On the Strength of Transverse 

Waves and Geometrical Detonation 
Cell Model for Gas Phase Detonations 

The MISTRESS User Manual H. H. Hilton 

The Chemical Shock Tube - R. A. Strehlow 

Implications of Flow Non- R. L^ Belford 

Idealities 


A. R. Zak 
C. D. Engel 

R. A. Strehlow 

A. Klavins 

R. A. Strehlow 
H. 0« Barthel 

J. R. Biller 


Journal Publication 

Progress in High Temperature 
' Physics "and Chemistry , 
Pergamon Press, N.Y., 2; 
127-176 (1968). 

J. Spacecraft, 5: 265-269 (1968) 


J. Spacecraft , 270-275 (1968). 
AIAA Journal, 7: 492-496 (1969). 


Progress in High Temperature 
Physics and CHCTrstiy 
Pergamon Press » N.y., 2: 

1-146 (1969). 

A, Klavins and L. H. Sentman 

Rev, Sci. Instr., 41: 1560-1567 (1970). 


R. A. Strehlow and J. R. Biller 
Comb, and Flame, 13; 577-582, 
' (1970). 
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Wave Structure 
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Investigation of a Two- 
Dimensionally Stiffened Panel 

AAE 69-5 On the Kinetic Equations for a 

Dilute, Short Range Gas 


AAE 69-6 The Sawtooth Column of the 

Supersonic Electric Arc in 
Sulfur Hexafluoride 

AAE 69-7 Theoretical and Experimental 

Analysis of Stiffened Panels 

Under Dynamic Conditions 

AAE 70-1 On the Interaction Between 

Chemical Kinetics and Gas- 
Dynamics in the Flow Behind a 
Cylindrical Detonation Front 

AAE 70-2 Preliminary Studies on the . 

Engineering Applications of 

Finite Difference Solutions of 
the Navier-Stokes Equations 

AAE 70-3 Some Aspects of the Surface Boundary 

Condition in Kinetic Theory 



A. Klavins 
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UILU-ENG 71 0503 


AAE 72-1 

UILU-ENG 72 0501 


Title 


Author 


A Study of the Transient Behavior of Fuel H. Krier 

Droplets during Combustion: Theoretical 

Considerations for Aerodynamic Stripping 

On the Solid Body Model for an Accelerating F. Klett 

Electric Arc 


Detonative Mach Stems 


R. A. Strehlow 
H, 0. Barthel 


An Investigation of Transient Phenomena R. J. Stiles 

in Detonations 


On the Role of Tangential Velocity Changes 
in the Scattering of a Molecular Beam from 
A Solid Surface 


C. C. Chrisinan 
L. H. Sentman 


Unconfined Vapor Cloud- Explosions - An R. A. Strehlow 

Overview 


Journal 

Publication 


with R. A. Strehlow, 

A. A. Adamczyk, 
Astronautica Acta , 

IJ;. 509-527 (1972) 

Chemical Physics 
Letters, 26:407-413(1974 


Fourteenth Symposium 
on Combustion, 1189- 
1200 (1973). 


AAE 72-2 
UILU-ENG 72 0502 


AAE 72-3 
UILU-ENG 72 0503 


Application of Illiac IV Computer to 
Numerical Solutions of Structural Problems 


On the Measurement of Energy Release Rates 
In Vapor Cloud Explosions 


H. H. Hilton 
A. R. Zak 
J. J. Kessler 
P. C. Rockenbach 

R. A. Strehlow 
L. D. Savage 
G. M. Vance 


Combustion Science 
and Technology , 
307-312 (1972). 


§ § 
II 
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UILU-ENG 

72 

AAE 73-1 


UILU-HNG 

73 

AAE 73-2 


UILU-ENG 

73 

AAE 73-3 


UILU-ENG 

73 

AAE 73-4 


UILU-ENG 

73 

AAE 73-5 


UILU-ENG 

73 

AAE 73-6 


UILU-ENG 

73 

AAE 73-7 


UILU-ENG 

73 

AAE 73-8 


UILU-ENG 

73 

AAE 73-9 


UILU-ENG 

73 





Journal 
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Author 

Publication 


A Performance Comparison o£ Several Numerical 

J. E. Pressing 


0504 

Minimization Algorithms 




Stresses and Damping in the Matrix of a 

A. R. Zak 


0501 

Composite Material 




Early Burning Anomalies in the XM 645 

U. Krier 

BRL Rept. No. 104 

0502 

Flechctto Cartridge 

D. R. Hall 

(1973) . 


Equivalent Explosive Yield of the Explosion 

R. A. Strehlow 


0503 

in the Alton Southern Gateway Yard, East St, 




Louis, 111., January 22, 1972 




Failure Studies of Gaseous Detonations 

R. J. Salm 

Acta Astronautica 

0504 



(in press) . 


An Investigation of Hydrogen-Oxygen-Argon 

J. R. Biller 


0505 

Detonations 




Interior Ballistic Predictions Using Data 

H. Krier 

Proc. 11th JANNAF 

0506 

From Closed and Variable- Volume Simulators 

S. A. Shimpi 

Combustion Meeting, CPI. 



M, J. Adams 

Publ. 261:17-30(1974). 


Theory of Rotationally Symmetric Laminar 

G. M. Vance 


0507 

Premixed Flames 

H. Krier 

Comb, and Flame J., 




365-375 (1974). 


Burning of Fuel Droplets at Elevated 

J. H. Rush 


0508 

Pressures 

H. Krier 

Comb, and Flame J. , 




377-382 (1974). 


An Impact Ignition Model for Solid 

H. Krier 


0509 

Propellants 

H. H, Hilton 

BRL Rept. No. 1707 



0. Olorunsola 

(1974). 



D, L. Reuss 
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ME 74-3 
UILU-ENG 74 

ME 74-4 
UILU-ENG 74 


ME 74-5 
UILU-ENG 74 

ME 74-6 
UILU-ENG 74 

ME 74-7 
UILU-ENG 74 
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UILU-ENG 74 
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UILU-ENG 74 
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Optimal Multiple- Impulse Direct Ascent 
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J. E, Trussing 
L. 11. Gross 

AIM J. , 12, 
885-889 (1974). 

0511 

The Structure and Stability of 
Detonation Waves 

R. A. Strehlow 


0501 

Model of Flame Spreading and Combustion 
Through Packed Beds of Propellant Grains 

H. Krier 

W, F. Van Tassell 
S . Raj an 
J. T. Ver Shaw 

BRL Report No. 147 
(1974). Int. J. Heat-Mass 
Transfer, 1377-86 (1975). 

0502 

0503 

On the Nature of Non- Ideal Blast 
Waves 

Viscous Incompressible Flow in Spiral 
Channels 

R. A. Strehlow 
A. A. Adamezyk 

W. F, VanTassell 

WSS/CI Paper 

No, 74-12, Pullman, Wash. 
(1974). 

0504 

Frequency Response Functions of a 
Disordered Periodic Beam 

J. N. Yang 
Y. K. Lin 

J. Sound and Vibration 
317-340 (1975) . 

0505 

Predicting Uniform Gun Interior Ballistics: 
Part I. An Analysis of Closed Bomb Testing 

H. Krier 
S. A. Shirapi 

Comb, and Flame J. 
229-240 (1975). 

0506 

0507 

Predicting Uniform Gun Interior Ballistics; 
Part II. The Interior Ballistic Code 

Predicting Uniform Gun Interior Ballistics: 
Part III. The Concept and Design of the 
Dynagun Ballistic Simulator 

H. Krier 
M. J. Adams 

H. Krier 
J. W. Black 

Proc. 11th JANNAF Comb. 
Meeting, CPIA Publ, 251: 
17-30 (1974). 

Proc. 11th JANNAF Comb. 
Meeting, CPIA Publ. 261: 
31-43(1974) . 

0508 

Process of Fluidization During Porous Solid 
Propellant Combustion 

H. Krier 
J. T. Ver Shaw 

AIM Paper 75-242 
(1975). 

0509 

An Analysis of Flame Propagation Through 
Coal Dust-Air Mixtures 

J. L. Krazinski 
H. Krier 

AIM Paper 74-1111 
(1974); 
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Sg 


An Interior Ballistics Prediction of the M549 

H. Krier 



0510 

Rocket' Assisted Projectile 

S. Shimpi 





E. Meister 


CJ ^ 


Dynamically Induced Thermal Stresses in 

A. Zak 


a 

0501 

Composite Material, Structural Panels 

W. Drysdale 


CO 


Numerical Analysis of Laminated, Orthotropic 

A, R. Zak 



0502 

Composite Structures 





The Characterization and Evaluation of Accidental 

R, A. Strehlow 

NASA CR 134779 

0503 

Explosions 

W. E. Baker 

G)une’ 1975 j. 

Also 




Progr. Energy ^ Comb 


Program Manual for the Eppler Airfoil Inversion 


Sc. (in press). 

0504 

Program 

W. G. Thomson 




Design of High Lift Airfoils with a Stratford 




0505 

Distribution by the Eppler Method 

W . G . Thomson 




Prediction of Flame Spreading and Pressure Wave 

H. Krier 

AIM J. 14: 


0506 

Propagation in Propellant Beds 


301-309 (1976) 


Vigorous Ignition of Granulated BedS' by Blast- 

H, Krier 

Int, J. Heat- 

•Mass 

0507 

Intact 

S. Gokhale 

Transfer 19: 





915-923 (1976) 


Solid Propellant Burning Evaluation .wit^ the 

H, Krier 

J. Ballistics 

1: 

0508 

Dynagun Ballistic Simulator 

T. G, Nietzke 
M. J, Adams 

103-149 (1976) 



J, W. Black 





E. E, Meister 




Structural Reliability § Minimum Weight Analysis 

H. H. Hilton 

.AlAA J, 


0509 

for Combined Random Loads § Strengths 


(in press) 
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Linear Viscoelastic Analysis with Random Material 
Properties 

H. H. Hilton 
J. Hsu 
d. S. Kirby 


0501 

0502 

Two Degree o£ Freedom Flutter o£ Linear Visco- 
elastic Wings in Two Dimensional Flow 

An Error Analysis o£ Computerized Aircra£t 
Synthesis 

C. F. Vail 
-H. H. Hilton 

V, V, Volodin 
H. H. Hilton 

In press AIAA J. 

In press 
J, of Aircraft 

0503 

Reactive Two-Phase Flow Models Applied to the 
Prediction o£ Detonation Transition in Granu- 
lated Propellant 

H. Krier 
M. Dimitstein 
S. S. Gokhale 

AIAA J. 16: 
177-183 (1978) 

0504 

Transient Temperature Response o£ Charring 
Composite Slabs 

J. E. Trussing 
H. Krier 

Int’l J, Heat 
Mass Transfer. 

519-522 (1978) 

0505 

0506 

Nonlinear Response o£ Laminated Composite 
Material Cylindrical Shells 

An Investigation o£ Blast Waves 
Generated from Non-Ideal Energy Sources 

A. R. Zak 
J. N. Craddock 

. - . A, -A. -Adamezyk 


0501 

Nonlinear Dynamic Analysis of Flat Laminated 
Plates by the Finite Element Method 

A, R, Zak 


0502 

An Investigation o£ Blast Waves Generated by 
Constant Velocity Flames 

R, T. Luckritz 


0503 

On the Blast Waves Produced by Constant Velocity 
Combustion Waves 

R. A. Strehlow 
R, D. Luckritz 


0504 

Direct Initiation of Detonation by Non-Ideal 
Blast Waves 

R, J, Cesarone 


0505 

The Blast Wave Generated by Constant Velocity 
Flames 

S, A. Shimpi 
R. A. Strehlow 
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AAE 77-15 
UILU-ENG 77 0515 
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Journal 

Author Publication 


Exploratory Studies of Flame and Explosion Quenching 


R. A. Strehlow 
L. C. Sorenson 
L. D. Savage 
H. Krier 


The Trajectory of a Liquid Droplet Injected Into the 
Wake of an Aircraft in Ground Effect 

Comparison of Viscoelastic and Structural Damping in 
Flutter 

The Blast Wave Generated by Constant Velocity Flames 


Wind Energy: History, Economics, and the Vertical 

Wind Turbine 

Final Report: Low Speed Airfoil Study 

Final Report: Propeller Study, Part I, Introduction 

and Overview 

Final Report: Propeller Study, Part II, The Design of 

Propellers for Minimum Noise 

Final Report: Propeller Study, Part III, Experimental 

Determination of Thrust § Torque on the Y0-3A Aircraft 


Direct Initiation of Detonation 


M, B. Bragg 
H. H. Hilton 

R. A. Strehlow 

R. T. Luckritz 
A. A. Adame zyk 

S. Shimpi 

T. R. Richards 


A. I. Orrasbee 

A. I, Orrasbee 

p o 

C. J. Woan 

ss 

S. A. Siddiqi 


K. R. Sivier 

A. I. Ormsbee 

R. A, Strehlow 
H. 0. Barthel 
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AAE 77-16 
UILU-ENG 77 0516 

The Effects of Energy Distribution Rates and Density 
Distribution on Blast Wave Structure 

R. A. Strehlow 
L. H. Sentman 


AAE 77-17 
UILU ENG 77 0517 

Modeling of Convective Mode Combustion Through Granu- 
lated Propellant to Predict Transition to Detonation 

H. Krier 
J. A. Kezerle 

17th Combustion 
Symposium: In Press 

AAE 78-1 

UILU ENG 78 0501 

Unsteady Internal Boundary Layer Flows with Applica- 
tion to Gun Barrel Heat Transfer and Erosion 

M. J. Adams 
H. Krier 


AAF 78-2 

UILU ENG 78 0502 

Extracting Burning Rates for Mul tip erf orated Propellant 
from Closed Bomb Testing 

H, Krier 


AAE 78-3 

UILU ENG 78 0503 

Lean Limit Flammability Study of Methane-Air Mixtures 
in a Square Flammability Tube 

J. Jarosinski 
R. A. Strehlow 


AAF 78-4 

UILU ENG 78 0504 

Interim Technical Report AFOSR 77-3336: ”An 

Investigation of the Ignition Delay Times For 
Propylene Oxide-Oxygen-Nitrogen Mixtures" 

E. E. Meister 


AAE 78-5 

UILU ENG 78 0505 

Final Report: A Distribution Model for the Aerial 

Application of Granular Agricultural Particles 

S, T. Fernandes 
A. I. Ormsbee 




